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Abstract: The structural and electronic properties of group V transition metal 
carbides: VC, NbC and TaC are studied using density functional theory with 
generalized gradient approximation for the exchange and correlation potential. 
Lattice constants, bulk moduli, elastic constants, energy band diagrams, density 
of states of the carbides are reported and their trends are discussed. From the 
band diagrams, the band separation, zone-centre d band splitting, non metal 
p and s band splitting, width of 2p band, and 4d band explains the insight of 
electronic structure of these compounds. 

Keywords: Density functional theory; energy band diagram; elastic 
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1. INTRODUCTION 


The transition-metal carbides belong to the family of refractory (hard 
metals). The refractory group can be formed by combining boron, carbon, 
or nitrogen with group-IV, V & VI transition series elements [13]. They 
are well known for possessing a number of extreme properties that are 
of scientifically interest and technologically important. Several of the 
compounds have high super-conducting transition temperature T c showing 
a strong electron phonon interaction [16], In the view of these properties, a 
large number of theoretical investigations have been carried out in the last 
few decades [1, 2, 7, 8, 12, 20, 21], The physical properties like extreme 
hardness are due to the presence of their unusual combination of bonding 
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mechanisms. The other properties like high melting points and electrical 
conductivity are of equal interest. These compounds often crystallize in 
rock salt structure [13]. The aim of this paper is to present the results of a 
theoretical investigations on the electronic structure and elastic properties 
of VC, NbC and TaC. The present work is organized as follows. The theory 
and method of electronic structure calculation are discussed in section 2. 
We present our results and discussion in section 3. In section 4, we draw 
some conclusions based on our results. 


2. METHODS OF CALCULATION 

The quantum many body problem with Bom-Oppenheimer approximation 
is much simpler, but still too far difficult to solve. Several methods exist to 
reduce this many particle Hamiltonian to an approximate one but in tractable 
form. One of the most powerful methods to reduce in such approximate form 
is Density Functional Theory (DFT) [19, 22, 25], Within this formalism, the 
total energy can be expressed as a functional of density of electron system, p 
and can be written by the following functional, 


E [p] = T 0 [p] + V H [p] + V xc [p] + V ext [p] 


where, T o [p] is the kinetic energy of a non-interacting electron system, V H [p] 
and V xc [p] are the Hartree and exchange-correlation contributions to the 
energy and V ext [p] is the energy due to the external potential of the system. 
Thus the corresponding Hamiltonian called the Kohn-Sham Hamiltonian is 


HKS=T 0 +V H +V xc +V ext 


2rm 


-Vr + 


47T£, 


■/ 


^f' + v xc + v„ 


The exact density of N electron system can be expressed in ground state by, 
N 

p( r) = E<A(T>* 0j(r) where the single particle wave functions 4>.(r) are the N 

lowest energy solutions of the Kohn Sham equation: H 4>. = Cdr of the N electron 
system. These equations can be solved self consistently in an iterative process. 
Thus the energy band is calculated using Full Potential-Linearized Augmented 
Plane wave (FP-LAPW) method within the Density Functional formalism in its 
local density approximation [17]. In this method, the unit cell is divided into non¬ 
overlapping atomic spheres (muffin tin spheres) and an interstitial region. Within 
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the muffin tin sphere, one constmcts the basis set as linear combination of radial 
functions multiplied by spherical harmonics. Whereas, a plane wave expansion is 
used as basis set in the interstitial region. The exchange and correlation potential is 
incorporated by using generalized gradient approximation (GGA) parameterized 
by Perdew et al [11], For every case the wave functions inside the MT spheres 
which are expanded into spherical harmonics are upto 1 = 9 and the RKmax = 10. 
The number of k points used for the integration procedure is 5000 which reduces 
to 165 irreducible k points inside the Brillion zone. In this work, the states treated 
as bands are: for V(3s,3p,3d,4s); Nb(4s,4p,4d,5s); Ta (5s, 5p, 5d, 6s); C(2s, 2p) 
and N(2s, 2p). The calculations are performed at the equilibrium lattice constants 
which are determined from a plot of total energy against the unit cell volume by 
fitting to the Mumaghan equation of state [10]. The overall simulation work is 
carried out using wien2k code [ 18]. 

3. RESULTS AND DISCUSSION 
3.1 Structural and Elastic Properties 

A series of different lattice constants are used to calculate the total energy and the 
corresponding unit cell volume to optimize the rock salt structure of the group V 
transition metal carbides (VC, NbC, TC). Figure l(a,b,c) shows the total energy 
curve as a function of unit cell volume and fitted with Mumaghan equation of 
state [10]. The obtained equilibrium lattice parameter (a Q ), bulk modulus (B 0 ), and 
pressure derivative of the bulk modulus (B’) are presented in table 1. These values 
are compared with previous theoretical and experimental results as calculated 
by other groups [3, 4, 6, 9, 13, 24], It shows that our results are reasonably in 
agreement with previous reported data. Hence this optimized structure is used 
for further study of elastic properties of the optimized lattice constant, the elastic 
constants C n , C p and C 44 are calculated using the stress-strain method [3,6,9]. 


Table 1: Lattice constant (a Q ), bulk modulus (B), pressure derivatives (B’) 
and elastic constants (C n , C 12 , C 44 ). 



VC 

NbC 

TaC 

a 0 (A°) 

4.165 (4.170) a 

4.683(4.470)° 

4.484(4.456) d 

B (Mbar) 

314.860(4.142) b 

300.031(302.000)° 

325.56(3.420)° 

B’ 

3.388 

3.894 

8.244 

C,,(GPa) 

5.056 

6.022 

8.129 

C p (GPa) 

1.288 

0.945 

0.819 

CJGPa) 

1.842 

1.743 

1.769 

a Ref.[25], b Ref.[17], 

c Ref.[ 13], d Ref.[ 11,10],' Ref.[l8]; 
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Figure 1: Energy verses Volume for (a)VC, (b) NbC, (c) TaC. 

The elastic constants C n , C p and C 44 are listed in table 1 and they satisfy 
the mechanical stability conditions for a cubic crystal system given by: 
(C n +2C p )> 0; C n C 12 > 0; C 44 > 0. Thus the group V transition metal carbides 
in rock salt structure is mechanically stable. 

3.2 Electronic Structure 

The energy band diagram of VC, NbC and TaC are shown in figure 2(a,b,c) 
respectively. In the energy band diagrams, one can see the lowest valence and 
conduction band. Also there is crossover of valence band and conduction band 
towards the Fermi level. It indicates that these carbides show metallic nature. 
Figure 3(a, b), Figure 4(a, b) and Figure 5 (a, b) show the partial density of 
states of V, Nb, Tc and C atom in VC, NbC, TaC respectively. To understand 
the contribution of orbital states in the band diagrams, we need to analyze 
thoroughly the partial density of states of the atoms in the system. The partial 
DOS, as shown in figure 3, figure 4 and figure 5 clearly shows that the valence 



Figure 2: Band diagram of (a) VC, (b) NbC , (c) TaC. 
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Figure 3: Partial DOS of (a) V in VC, (b) C in VC. 




Figure 4: Partial DOS of (a) Nb in NbC, (b) C in NbC. 




Figure 5 : Partial DOS of (a) Ta in TaC, (b) C in TaC. 

band in energy band diagrams are mainly contributed by the nonmetal 2s 
orbital state of C atom. The valence bands are mainly due to nonmetal p state 
and metal p state. 
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Table 2: Band separations (Ry) in 
LDA calculations. 

band diagram obtained from GGA, APW, 


VC 

GGA a 

APW b 

NbC 

GGA a 

GGA C 

APW C 

LDA b 

TaC 

GGA a 

LDA d 

Thapa, RK 

e f 

0.876 

—- 

1.020 

1.055 

0.770 

1.085 

1.143 

0.091 



0.199 

0.210 

0.174 

0.173 

0.250 

0.138 

0.129 

0.081 



0.152 

0.070 

0.110 

0.119 

0.133 

0.128 

0.082 

0.084 


E d - E P 

0.001 

0.100 

0.030 

0.031 

0.150 

0.032 

0.093 

0.356 


E -E 

1.031 

0.940 

0.965 

0.965 

0.870 

0.969 

0.981 



P s 


a Ref. Our result, b Ref. [5], c Ref. [14 ,23], d Ref. [15] 

It is also notice that there is a small contribution of metal d states. One can 
observe that the 2s band has maximum dispersion between T and L. Above 
this band, there lies the 2p non-metal and 4d metal states. These bands overlap 
and mixed in A. The 4d metal state is further decomposed to t, and originating 
from T 25 . and c. originating from T 12 . So the width of the d band is given by 
AE d = E(T 12 ) - E(T 25 ,). In our calculation, the energy gap E gL is the band 
separation of lowest band and valence band at L point. The zone centre d band 
splitting AE d , the metal-d & non-metal-p splitting, E d - E = [ 2/5 E (T 12 ) + 3/5 
E(T 25 ,)] — E(T 15 ), nonmetal p & metal s splitting, E p - E s = E (T 1S ) — E(Tj) 
are given in table 2. Our band structure calculation parameters are in good 
agreement with the reported in the literatures [2,5, 14, 15, 23], 

CONCLUSIONS 

This work reports the results of a systematic study on the electronic and 
structural properties of group V transition metal carbides. The optimized 
lattice constant of VC, NbC and TaC are found as 4.165A°, 4.683 A° and 
4.4834 A° which are in good agreement with the reported experimental and 
calculated values. The elastic constants (C u , C , C 44 ) satisfy the mechanical 
stability conditions: (C n +2C 12 )> 0; C n C 12 > 0; C 44 > 0. Thus the group V 
transition metal carbides (VC, NbC, TaC ) are found to be stable compound in 
rock salt structure. From a comparison of the band structures, one concludes 
the behavior is directly related to the hybridization of the metal-d and non 
metal-p valence states near E . The change in the overlapping and mixing 
between metal d (3d for V, 4d for Nb and 5d for Ta) and non metal-2p bands 
are random. The value of band gap, zone-centre d band splitting, non metal 
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p and s band splitting explains the insight of electronic structure of group V 
transition metal carbides. 
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